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ABSTRACT 

Exploiting a five-year span of data, we present improved timing solutions for the five millisecond pulsars 
known in the globular cluster NGC 6752. They include proper motion determinations for the two outermost 
pulsars in the cluster, PSR J1910-5959A and PSR J1910-5959C. The values of the proper motions are in agree- 
ment with each other within current uncertainties, but do not match (at Act and 2a level respectively) with the 
value of the proper motion of the entire globular cluster derived in the optical band. Implications of these re- 
sults for the cluster membership of the two pulsars are investigated. Prospects for the detection of the Shapiro 
delay in the binary system J1910-5959A are also discussed. 

Subject headings: pulsars: individual (PSR J1910-5959A; PSR J1910-5959B; J1910-5959C; J1910-5959D; 
J1910-5959E) — globular clusters: individual (NGC 6752) 



1. INTRODUCTION 

The globular cluster (GC) NGC 6 752 is known t o host 
five millisecond pu lsars (MSPs) JD'Amico et alJ 1200 It 
iD'Amico etaTJl2002l hereafter Paper I). PSR J1910-5959B 
and PSRJ1910-5959E (hereafter PSRB and E respectively) 
reside in the central region of the GC and show large neg- 
ative P values, which are interpreted as an effect of the GC 
gravitational potential well (Paper I). This in turn implies a 
large mass-to-light ra tio in the central region of NGC 6752. 
iFerraro et alJ J2003bl) recalculated the center of gravity and 
studied the luminosity profile of this cluster: combining 
their HST data with the P value of PSRB and PSRE, they 
put a firm lower limit on the central mass-to-light ratio of 
M/L v > 5.5M Q /L . Also PSR J1910-5959D (PSRD) is 
located close to the GC center. Its P value is positive and 
of the same order of magnitude of PSR B and E, suggesting 
that also for PSR D the P value is dominated by the gravita- 
tional potential well (Paper I). PSR J 1 9 1 0-5959C 7 (PSR C) is 
located at a projected distance 9± = 2'.6 from the GC cen- 
ter (Paper I), which i s much larger than th e cluster's core ra- 
dius r c = 5'.'2±2'.'4 (Ferraro et al. 2003b). The binary pulsar 
J1910-5959A (PSR A) is located at an even larger distance 
from the GC center (9± = 6'A, Paper I), the largest offset 
known for a GC pulsar. 

The positions of PSR A and PSR C are unexpected since 
mass segregation should have driven the two neutron stars 
close to the GC center in a time scale (< 1 Gyr) much shorter 
tha n the time s i nce their for mation (~ 10 Gyr). In particu- 
lar, iColni et alJ 120021 12001 explored various scenarios to 
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explain the unusual position of PSR A, invoking a dynami- 
cal encounter in the inner region of the GC. The most prob- 
able picture is that PSR A was originally in the GC central 
regions and it has been expelled to the outskirts by the in- 
teraction with either a single massive black hole (BH) or a 
binary [BH + BH] of unequal mass. Timing results in Paperl 
indicated a low-mass white dwarf as the most probable com- 
panion for PSR A. This has bee n confirmed by B assa et in 
(2003) and Ferr aro et al1J2003al) . who identified with Hubble 
Space Telescope observations the companion of PSR A with a 
helium white dwarf star whose mass is M co ~ 0.17-0. 20 M 
and whose photometric properties are compatible with its be- 
longing to NGC 6752. 

The issue of the association of PSR A to NGC 6752 has 
been recently revisited using spectroscopic observations of 
the op tical companion to th e pulsar, performed with the ESO- 
VLT. [Cocozza et al. (2006) found full agreement (at la) be- 
tween the radial velocity of the center of mass of the binary 
7 = -28.1 ±4.9 km s -1 and the overall c luster radial v elocity 
v'6752 = -27.9 ± 0.8 km s" 1 , obtained bv lHarrisI i 19961) (cata- 
log revision 2003) averaging various determinations. This is a 
strong indication in favour of the association of the pulsar with 
NGC 6 752. However, using the same data set, Bas sa et all 
(2006) compared the systemic velocity of the binary with that 
of nearby stars which certainly belong to the cluster and con- 
cluded that they are only marginally consistent at 2cr level. 

In this paper we present timing results based on more than 
five years of regular observations. In particular, with a much 
longer available data span we have been able to measure 
proper motions of PSR A and PSR C. The new timing solu- 
tions as well as the pulse profiles for all the millisecond pul- 
sars are presented in § 2, § 3 reports on the proper motion de- 
terminations and implications for the cluster membership of 
the two pulsars are discussed in § 4. 

2. OBSERVATIONS AND IMPROVED TIMING PARAMETERS 

Regular pulsar timing observations of NGC 6752 have been 
carried out since September 2000 with the Parkes 64 m radio 
telescope at a central frequency of 1390 MHz, using the cen- 
tral beam of the multibeam receiver or the H-OH receiver. 
The hardware sy stem is the same as th at used in the discov- 
ery observations (D'Amico et al. 2001). The effects of inter- 
stellar dispersion are minimised by using a filterbank having 
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5 12x0.5 MHz frequency channels for each polarization. Af- 
ter detection, the signals from individual channels are added 
in polarization pairs, integrated, 1 bit-digitized every 125/xs 
(80 /is in recent observations), and recorded to magnetic tape 
for off-line analysis. Pulse times of arrival (TOAs) are de- 
termined by fitting a template profile to the observed mean 
pulse profiles and analysed using the program TEMPO 8 and 
the DE405 solar system ephemeris. 

Tabled summarizes the best fit values and uncertainties 
(chosen to be twice the nominal TEMPO errors) for the pa- 
rameters entering our timing solutions, whose residuals are 
displayed in Figure^ The same figure presents a high signal- 
to-noise profile obtained for each of the pulsars by folding the 
best available data according to the reported ephemerides. 

The new positional and rotational parameters at the refer- 
ence epoch are all compatibile with those reported in Paper I 
(assuming 3er uncertainties for the values quoted in Paper I). 
However, the MJD range of the available TOAs is now ^3.5 
times longer than for Paper I and hence the accuracy of the so- 
lutions has improved correspondingly. Orbital parameters for 
PSR A, obtained using the ELLl model of TEMPO, have also 
been measured with a higher precision than in Paper! Figure 
|2]shows that no trend is evident in the timing residuals plotted 
with respect to the orbital phase for the timing solution given 
in Table [2 An additional constraint on the orbit of PSR A 
has resulted from the recent optical observations of the pulsar 
comp anion. Spectroscopy dCocozza et alJl2006t iBassa et alJ 
2006) has provided us with a measurement of the mass ratio 
in the binary, whereas multi-color photometry h as set the pos- 
sible range for the mass of the secondary star tFerraro et alJ 
I2003at IBassa et al] l2006). Combining these results gives a 
limit on the orbit al inclination > 70° dCocozza et alJl2006t 
Bass a et all2006T) . 

The size of the expected Shapiro delay is nominally larger 
than the rms residual of the timing solution (see Table Q for 
any i > 70°, but, except for inclination angles near 90°, a large 
part of the Shapiro delay is absorbed in the Roehmer delay 
( Lang e et alJ2001l) . In fact, no clear trend is visible in the tim- 
ing residuals even after binning the TOAs in orbital phase (see 
Figure |2j, indicating that the magnitude of the unabsorbed 
component of Shapiro delay is below the present uncertainty 
in the TOAs. Therefore it is not surprising that fitting the 
available TOAs with TEMPO has not led to any significant de- 
termination of the Shapiro parameter s. Inspection of Figure 
|2]also shows that the present uncertainties on the TOAs allow 
us only to exclude very extreme orbital inclinations i > 89°. 
Simulations show that a factor ~ 2- 3 improvement in timing 
precision is needed in order to obtain a useful constraint on s. 
This will require an additional ~ 10 years of observation with 
the present instrumentation and collection rate of TOAs. 

The still unassessed effects of Shapiro delay may also af- 
fect the new determination of the binary eccentricity, for 
which in Paper I only an upper limit was available. Neglect- 
ing Shapiro delay, the measured value is e = 3.4(12) x 10~ 6 
(here and everywhere in this paper the errors are quoted at 
twice the nominal rms values given by TEMPO). However, 
for 70 < i < 89° and 0.17 M < M co < 0.20 M an unmod- 
eled Shapiro delay can introduce an apparent eccentricity in 
the range 1 - 3 x 10~ 6 . The determination of e must still be 
considered provisional and e — 4.6 x 10~ 6 is a reliable upper 
limit. 

8 see http://www.atnf.csiro.au/research/pulsar/timing/tempo 



The small eccentricity of PSR A's binary system is typical 
of fully recycled binary millisecond pulsars and is consistent 
with the ef fects of ran dom encounters with other cluster stars 
(Rasio & Hegsie 1995). The upper limit on e is also compat- 
ible with the offset position of PSR A resulting from an inter- 
action which occurred ~ 1 Gyr ago between the already re- 
cycled binary system including PSR A and a WD companion 
[PSR A +WD1 with a bina ry black hole of a few tens of solar 
masses (Colpi et al. 2003, see § 1). We note that the value of 
e also fits with the hypothesis (Ba ssa et"al]l2003t IColpi etafl 
120031) that a dynamical encounter with a single BH, whose 
mass is higher than a few hundred M Q , may have simultane- 
ously ejected the progenitor of [PSR A+WD] and triggered 
the recycling process in the binary, which in turn circularized 
the system and removed any information on its post-encounter 
eccentricity. However, the value of e does not agree with an 
ejection event involving the already formed [PSR A + WD] bi- 
nary and a single BH. In this case, Colpi et alJ (1200 3) showed 
that the post-encounter eccentricity of [PSR A+WD] would 
be significantly larger, up to values of 10~ 4 - 10~ 2 and only 
slightly affected by subsequen t random encou nters with nor- 
mal stars of the cluster (Rasio & Heasie 1995J). 

The mean flux densities at 1400 MHz (Shoo) in Tableware 
average values, derived from the system sensitivity, the ob- 
served signal-to-noise ratio, the shape of the pulse profile, 
the displacement of the pulsars with respect to the center of 
the telescope beam and assuming flux density values corre- 
sponding to half the detection limit for the non-detections due 
to the strong interstellar scintillation effects on the pulsars in 
NGC 6752 (see § 3). The uncertainties on the values of S1400 
may reach ~ 30% f or the faintest sourc es. For a distance 
d = 4.45 ±0.15kpc ( iGratton et "all 120031) the inferred radio 
luminosities at 1400 MHz of the two millisecond pulsars in 
the cluster's outskirts are Z4400 = Snood 2 ~ 4-5mJykpc 2 , 
a value in the middle of the distribution of the luminosities of 
the millisecond pulsars in 47 Tucanae (Camilo et al. 2000). 

3. PROPER MOTION DETERMINATIONS 

The main improvement in our timing solutions is that 
proper motion determinations for the two outermost pulsars 
in NGC 6752 are now available. In Table ^proper motion 
components in right ascension and declination are reported as 
well as the corresponding proper motion amplitude and po- 
sition angle (PA, measured counterclockwise from north to- 
ward east). Proper motion uncertainties depend on the length 
of the data span and on the number, the degree of uniformity 
and the errors of the TOAs along the data span. The differ- 
ent precisions in our measurements are mainly due to the dif- 
ferent number of high quality TOAs available for each pul- 
sar, as shown in Table \l\ Measurement of good TOAs for 
the faintest pulsars is possible only when interstellar scintilla- 
tion enhances their signal: this is the reason that in the timing 
analysis of PSRs B, D, and E we used a significantly smaller 
number of TOAs than for PSR A. The flux density of PSRC 
is similar to that of PSR A and the effects of interstellar scin- 
tillation are also comparable. The difference in rms residual 
between the timing solutions for these two pulsars is primar- 
ily due to the different pulse widths, being ~ 7 times larger 
(at 50% of the peak) for PSR C than for PSR A. 

Figure|3] presents a geometrical representation of the ex- 
pected motion in the plane of the sky (during next 10 4 yrs) 
of PSR A, PSRC, and of the center of NGC 6752, as de- 
rived from their measured proper motions. The proper mo- 
tion for the center of the globular cluster has been obtained 
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bv lDinescu et alJ i 19991) . by comparing two optical observa- 
tions taken 25 years apart. The values for the components are 
H a cosS = -0.7 ±0.8 mas yr -1 and fig = -2.9±0.9masyr _1 . 
Their derivation required a transformation of the coordinate 
system at the epoch of the first observation to the coordinate 
system at the epoch of the second observation and the use of 
distant field galaxies as reference. The inset in Figure|3]shows 
a comparison between the proper motion vectors of PSR A 
and PSR C (with their uncertainties), and the optical proper 
motion vector of the cluster. The proper motions of PSR A 
and PSR C are compatible with each other, but they are not in 
agreement with the optical proper motion of NGC 6752, at 4a 
and 2a confidence levels respectively. 

4. DISCUSSION 

Since the escape velocity from a globular cluster is usu- 
ally significantly lower than the typical transverse velocity of 
these stellar systems, it is expected that the proper motion of 
a cluster pulsar will largely reflect the overall motion of the 
cluster. For NGC 675 2, the escape velo city from the central 
region is ~ 30km s -1 (iColpi et al.12 003) and the space veloc- 
ity is ~ 62 km s -1 with respect to the Solar Syst em barycenter, 
based on the proper motion measurement by Dinescu et al. 
( 1999) and the dista nce derived from the distance modulus 
(Gratt oiTet al.ll200 3f). Observations over a much longer data 
span may reveal the peculiar (orbital) motion of a pulsar in 
the cluster's gravitational potential well. 

Is it possible that the discrepancy between the proper mo- 
tions of PSR A and PSR C and the optical proper motion of 
NGC 6752 (§ 3) could be an indication that the two pulsars 
are not associated to the globular cluster? 

In Paper I it was estimated that the probability 9 for PSR A to 
be a Galactic field millisecond pulsar superposed by chance to 
NGC 6752 (at a distance of 6'A off its center) is of order 10" 5 . 
The compatibility of the measured proper motions of PSR A 
and PSR C reinforces the unlikeliness for these two MSP to be 
Galactic field objects by chance superimposed to the globular 
cluster. 

Assuming that both PSRs A and C are members of 
NGC 6752, the discrepancy between pulsar and GC proper 
motions, measured in the radio and optical band respectively, 
may result from the different methods used for determining 
the proper motions in the two spectral bands. In fact, simi- 
lar discrepa ncies have already been noted for the pulsars in 
47Tu canae (iFreire et alJl200lL 120031). in M4 (TThorsett et al.l 
1999) and, more recently, in Ml 5 dJacobv et all2006l) . 

However, the discrepancies in these clusters may not 
easily be ascribed to a common systematic effect affect- 
ing all the optical measurements. The optical proper mo- 
tion for 47Tucanae was directly measured based on Hip- 
parcos observations. The prop er motion determination for 
M4 (|Cudworth & Hanson! Il993l) was based on the determi- 
nation of its motion relative to a set of reference field stars 
whose proper motion relative to the Sun has been in turn ob- 
tained by combining their position in the Galaxy (through 
their parallax) to a dynamic model for the nearby regions 
of the Galaxy where these referen ce stars reside. Fin ally, 
the proper motion for NGC 6752 tDinescu et alJ U997T) re- 
sulted from the comparison of two photographic plates taken 

9 It is worth nothing that this probability does not account for the similar 
values of the dispersion measure of PSR A and PSR C. Given the uncertainty 
in the Galact ic electron layer, it is difficult to quantify the probability for this 
coincidence iBassaetal. 2006). However, it certainly further decreases the 
total probability for a chance superposition. 



25 years apart, using distant field galaxies as reference ob- 
jects. In the case of M15, four different optical determina- 
tions have been per formed (Geffert etai]fl993l: IScholz et alJ 
119961 lOdenkirchen et alJll997t ICudworth & Hansonlll993l) . 
three of which are incompatible with the pulsar proper mo- 
tions. These three non-matching measures were respec- 
tively based on the c omparison be tween photographic plates 
at different epochs JGeffert et alJ 119931 Ischolz et all 1 19961) 
and t he use of reference sta rs from Hipparcos observa- 
tions ( Odenkirchen et alJl!997l) . Only the measurement by 
Cudworth & Hanson ( 1993) is in agreement with the appar- 
ent motions of th e three pulsars investigated. However it is 
worth noting that Cudworth & Hanson! d!993l) measured the 
optical proper motion for M15 applying the same method as 
was used for M4, which in that case led to a discrepant proper 
motion. 

The discrepancy may be alternatively ascribed to very fast 
peculiar motions of PSR A and PSR C inside the cluster grav- 
itational potential well. Assuming that both the pulsar proper 
motions and the optical proper motion of the cluster are cor- 
rect, the relative 2-D velocity vectors of the pulsars with re- 
spect to the cluster center are roughly directed towards the 
cluster inner regions, as is indicated in Figure This would 
mean that [PSR A+WD] cannot be now in the phase of ejec- 
tion from the cluster and that it is not at its farthest distance 
from the GC center along its orbit insi de the cluster gravi- 
tational well. For d = 4.45±0.15kpc JGratton et alJl2003l) 
the relative transverse speed of PSR A would be V re i,A = 
51 ± 16kms _1 . NGC 6752 can provide a gravitational pull 
strong enough to retain PSR A at its actual location with a 
peculiar velocity V re i,A only if the mass enclosed within the 
pulsar projected position is M enc ] > 1.18 x 1O 6 M . This is 
in contrast to the total mass value for t he cluster obtained 
with HST observation (Sabbi et al. 2004), which is lower by 
an order of magnitu de. Using again the distance mod ulus in 
iCgatton et alJ (120031) . the apparent magnitude giv en bv [Harrisj 
( 1996) and the colour excess £b-v = 04 in Ferra roet alJ 
( 1999), the resulting overall mass-to-light ratio would be 
M enc i/L > 8.4M /L , which is unreasonably high for a GC, 
unless we admit an initial mass function much flatter than usu- 
ally estimated (so producing a very large number of under- 
luminous stellar remnants) or the presence in the cluster of 
a significant amount of dark matter. At a more conservative 
confidence level (4a) for the relative velocity V re i,A, the mass 
to light ratio would result M enc \/L > 2.5M Q /L Q , again im- 
plying a dynamic mass much higher than the mass derived 
from optical observations. 

A further test of the cluster membership of PSR A and 
PSRC will be possible in future. It will involve the compar- 
ison of the proper motions of PSR A and PSR C with those 
of the three pulsars close to the cluster core, whose asso- 
ciation with NGC 6752 is unambigously proved by the very 
strong gravitational pull affecting the values of their spin pe- 
riod derivative. This task will take some years: our simula- 
tions show that, with the present accuracy and collection rate 
of the TOAs and if the three innermost pulsars display the 
same proper motion as PSR A, a 3a determination will require 
a total data span of about 8-10 years. 
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FIG. 1 . — Fit residuals (left panels) and pulse profiles (right panels) for the five pulsars known in NGC 6752. Mean pulse profiles shown in the right panels are 
the sum of the observed profiles with the highest signal-to-noise ratio. The adopted binning (64 bins) matches the time resolution of the profiles. 
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FIG. 2. — faj Fit residuals versus orbital phase for PSR A obtained from the timing solution of Table fTI (b) Timing residuals binned in 42 orbital bins. The 
central values and the plotted uncertainties result from a weighted average (and error propagation) performed on all the available TOAs in each orbital bin. The 
lines represent the expected trends of the timing residuals when the Shapiro delay is not included in the timing model. The two upper curves are for an orbital 
inclination i = 89°, whereas the two lower curves are for i = 80°. The mass of the companion star is taken as 0.20 Mq (solid line) and 0.17 M (dotted line) in 
both sets of curves. 
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FIG. 3. — Main panel - Positions and expected changes (assuming uniform motion) after 10 years for PSR A, PSRC and the center of NGC 6752 relative to 
the present position of the cluster center. The uncertainties in the expected final positions are described by boxes whose size is given by the propagation of the 
uncertainties on the proper motions in right ascension and declinati on (2a confidence lev el). Proper motion uncertainties for the pulsars are from TablefTI while 
the uncertainties for the optical proper motion of the cluster are from Dinescu et al. 1 1999). The dashed circle represents the portion of the cluster enclosed within 
the half-mass radius i\ m = l'.9 i Trager et al. 1993). Lower-left panel - Comparison of the motions of the two outermost pulsars (PSR A: solid line; PSRC: dotted 
line)and the globular cluster (dashed line) relative to their present position. 
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TABLE 1 

Measured and derived parameters for the pulsars in NGC 6752. 



Parameter 


PSRA 


PSRB 


PSRC 


PSRD 


PSRE 


R.A. (J2000) 


19:11:42.75562(8) 


19:10:52.0556(5) 


19:11:05.5552(4) 


19:10:52.4163(5) 


19:10:52.1572(6) 


Decl. (J2000) 


-59:58:26.904(1) 


-59:59:00.861(6) 


-60:00:59.700(4) 


-59:59:05.479(5) 


-59:59:02.087(7) 


li a cos 8 (mas yr ) 


-3.3(2) 


- 


-4.1(17) 


- 


- 


fi s (masyr 1 ) 


-3.6(3) 


- 


-4.6(25) 


- 


- 


(mas yr -1 ) 


4.8(3) 


- 


6.2(22) 


- 


- 


PA fl (deg) 


222(3) 


- 


221(20) 


- 


- 


P(ms) 


3.2661865707908(1) 


8.357798500844(2) 


5.2773269323093(15) 


9.035285247765(4) 


4.571765939750(2) 


P b (ss-') 


2.947(2) xl0~ 21 


-7.9041(5)xl0- 19 


2. 16(2) xKr 21 


9.6431(6)xl0" 19 


-4.3435(3)xl0- 19 


Epoch (MJD) 


51920.0000 


52000.0000 


51910.0000 


51910.0000 


51910.0000 


DM (pc cm -3 ) 


33.705(3) 


33.33(6) 


33.29(5) 


33.28(2) 


33.31(3) 


P orb (days) 


0.8371134769(1) 


- 


- 


- 


- 


a sin i (1-s) 


1.2060461(8) 


- 


- 


- 


- 


T mc (MJD) 


51919.20647998(16) 


- 


- 


- 


- 


esmu) 


3.3(12) xl0~ 6 










ecosw 


0.9(13)xl(r 6 










f(M c ) (Ms) 


0.002687854(6) 










M £ .,„,„, C (M Q ) 


0.19 










MJD Range 


51710-53836 


51745-53769 


51710-53836 


51745-53731 


51744-53836 


Number of TOAs 


450 


44 


246 


124 


70 


r.m.s. residuals (fis) 


5.0 


18 


29 


24 


25 


Offset'' (arcmin) 


6.37 


0.06 


2.56 


0.05 


0.05 


Wio e @ 10% (ms) 


0.6 


1.3 


2.8 


1.1 


1.1 


W 50 f @50% (ms) 


0.4 


0.6 


1.3 


0.7 


0.6 


5i4oo (mJy) 


0.21 


0.05 


0.24 


0.05 


0.07 



"Position angle of the proper motion vector. 

b As discussed in Paper I, the observed P of PSRB, PSRC, PSRD, and PSRE are strongly affected by the gravitational potential well of the globular cluster. 
Useful constraints on the instrinsic spin-down rate Pj can hence be inferred only for PSRA. Correcting the observed value of P for ((') the Galactic differential 
rotati on and the vertical acceleration in the Galactic potential (see e.g. Damour & Tavlot 1991), for (ii) the centrif ugal acceleration of the pulsar i Shklovskii 
1970), and for (Hi) the contribution of the cluster potential well (estimated according to the recipe of Phinnev 1 1992) and using the luminosity density profile of 
NGC 6752 published by Ferrar o et alj I2003bl) ) gives P, < 6 x 10 21 s s 1 . We have also adopted M/Ly =5.5 i Ferraro et al. 2003b) in order to obtain the firmest 
upper limit on the intrinsic spin-down rate of PSR A. This translates to a lower limit to the pulsar spin-down age of 0.5(P/P,) ~ 8.6 Gyr and upper limits to the 
surface dipole magnetic field of 3.2 x 10" \J PP] ~ 1 .4 x 10 8 G and to the spin-down luminosity of 4tt 2 IP]/P 3 = 6.9 x 10 33 erg s _1 (/ being the moment of inertia 
of the neutron star, set equal to 10 45 g cm 2 ). 

c 'The minimum mass is calculated assuming a pulsar mass of 1 .35 M, ., and an inclination for the orbital plane with respect to the line of sight of 90° . 

d The offset of the pulsars is calculated with respect to the position of the cluster's center of gravity reported by Ferraro et al. 1 2003b). 

e Pulse width at 10% of the height of the main peak. 

f Pulse width at 50% of the height of the main peak. 



